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Over the last two decades, the HERA collider has provided a large 
amount of new information about QCD dynamics at high energy. While 
the most appreciated are the measurements of the proton structure func- 
tions in a wide range of parton momentum x and virtuality Q 2 , it is hard to 
believe that some of the observations at HERA which do not fit the simple 
picture of DGLAP dynamics would not get amplified at the LHC, possi- 
bly rendering certain approaches to searches beyond the Standard Model 
inadequate. 
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1. Introduction 

Progress in the next decade (and most probably even longer) in High 
Energy Physics will be totally overshadowed by the results expected from 
the Large Hadron Collider (LHC) at CERN (Geneva) where high energy 
protons will collide at an unprecedented center of mass energy of 14 TeV, 
any time soon. For quite some time, theorists claim that the success of the 
Standard Model (SM) of the unified electroweak interactions and of strong 
interactions cannot be understood, unless it is a low energy realization of a 
more sophisticated construct, which should manifest itself already at TeV 
scales. It still remains to be seen whether the Standard Model Higgs, whose 
mass is constrained by precision electroweak measurements pQ , is realized in 
nature. A low mass of the SM Higgs, about 120 GeV, creates the so-called 
hierarchy problem, that is the need at high energy for a tremendous fine- 
tuning of quantum-loop corrections to keep the Higgs mass light enough for 
it to be relevant for electroweak symmetry breaking. The favored candidates 
for the new horizon are supersymmetric (SUSY) theories and/or theories 
with universal extra-dimensions (UED). There are many realizations of these 
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theories (see [2] for a recent overview). Their common denominator is the 
appearance of new particles or states, which will be signaled by a departure 
of the measured cross sections from expectations of the SM. 

The generic LHC signatures of SUSY models are isolated leptons, few 
energetic jets and, for iZ-parity conserving SUSY, large missing transverse 
energy [3J. Similar signatures may be expected in the various realizations of 
the UED theories [HE]. These signatures may also originate from the SM, 
dominated by hard QCD interactions. In the context of LHC physics, QCD 
constitutes the background (though for many it is in itself a interesting 
laboratory for the dynamics of high energy QCD). The question is then 
how reliable is the background estimate. This is where the QCD studies 
at the HERA ep collider become very relevant, though the energy scales 
are mostly lower than the ones expected at LHC. The same could be said 
of the Tevatron pp collider studies, if not for the fact that the problems 
encountered at the Tevatron can only become worse at the LHC. And this is 
because of the composite nature of hadrons, nucleons in this particular case, 
which have a rich dynamical internal structure, that cannot be derived from 
perturbative QCD. This structure is directly probed in the deep inelastic 
scattering (DIS) of leptons on hadrons. 



2. Hard pp interactions 

The reactions that will be mostly investigated at the LHC are hard pp 
interactions in which partons in the proton, quarks and gluons, interact 
producing either high transverse momentum partons that materialize in 
the form of hadronic jets or new heavy particles which decay into leptonic 
and/or hadronic final states. The cross section for the hard interactions can 
be calculated in perturbative QCD, though to a limited precision due to 
higher order corrections. The hadronization process is soft in nature and 
requires modeling. Due to the complicated nature of resulting final states, 
even the expectations from pure QCD processes require MC simulation, 
which limit the precision with which the QCD background can be estimated. 

There are the following four components in modeling hadronic final 
states in hard pp collisions [B] : 

1. the hard subprocess at the parton level. At leading order it is a 
2 — > 2 process and the cross section is expressed through 

^? = E/^i /<*»*/ ^(xi,g 2 )/,(x 2 ,Q 2 )^g^5 [pi - f fj , 

(1) 

where pt is the transverse momentum of the outgoing partons k and I and 
s = x\X2S is the squared center of mass energy of the partonic scattering, 
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with s the pp center of mass energy squared and Xij the fraction of the 
proton momenta carried by the incoming partons i and j. The cross section 
for parton scattering is given by a, calculable in perturbative QCD, and 
the probability of finding the partons in the proton is given by the parton 
density functions (PDF) fi{x,Q 2 ), where the hard scale Q 2 is assumed to 
be Q 2 = p%. The variables t, u and s are the usual Mandelstam variables; 

2. initial and final state gluon radiation from the interacting partons, 
which is meant to mimic higher order QCD processes; 

3. additional semi-hard or hard interactions between the remaining 
partons. They arise naturally in that the integrated cross section given 
by (H|) diverges for low p? values and <7i n t becomes bigger than the total 
cross section, unless multi-parton interactions (MPI) are introduced; 

4. hadronization, the formation of hadrons from partons and from beam 
remnants. 

Each of these steps is subject to modeling and various Monte Carlo 
generators 13 QUI [11] adopt different approaches, in particular for 
the last three steps. Each one of the modeling steps is accompanied by 
uncertainties that are hard to quantify. In the following, the discussion 
will focus on the extraction of parton density functions which are essential 
for the calculation of the cross sections and on the validity of the DGLAP 
evolution in describing the QCD dynamics at high energy. 



3. Parton density functions in the proton 

There are many parameterizations of the proton PDFs derived in global 
fits of the NLO DGLAP evolution equations to the proton structure func- 
tions and hard scattering measurements. The most commonly used are the 
ones provided by the MRST(MSTW) [12] and the CTEQ groups p3]. They 
incorporate not only the DIS data, but also data from hadron-hadron inter- 
actions. A big effort is now vested in quantifying the uncertainties |14j for 
the precision measurements at the LHC. The worrisome part is that these 
uncertainties usually come out smaller than the differences between the 
various groups and this is particularly true with the gluon distribution, the 
contribution of which dominates the QCD cross section at LHC. Recently, 
the two HERA experiments, HI and ZEUS, provided the most precise PDFs 
by combining their HERAI measurements of the F2 structure function of 
the proton [15]. The "cross-calibration" reduces dramatically the system- 
atic uncertainties, as shown in Fig. [TJ Soon, one may expect even better 
precision, when the high luminosity measurements of neutral current and 
charged currents e^p interactions from HERAII become available. 
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HI and ZEUS Combined PDF Fit 



Q 2 = 10 GeV ! 




Fig. 1. Left - comparison of the HI and the ZEUS PDFs as a function of x at 
Q 2 = 10 GeV 2 and right - the result of fitting PDFs to the combined measurements 
of i<2 at the same Q 2 value. 



3. 1 . Low x regime 

The extraction of PDFs assumes the validity of the DGLAP evolution 
equation in the full range probed by the data, including the very low x 
regime of HERA. The only marker justifying this approach is the good x 2 
of the global fit. However many parameters are involved and the lever arm 
at low x is relatively small (see Fig. [2]) . A large fraction of the low x cross 
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Fig. 2. The coverage of the x and Q 2 plane by existing measurements and the 
expected coverage in the LHC pp experiments. 
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section in ep scattering is due to heavy flavor production, in particular cc 
pairs. A compilation of HERA measurements [16J of F£ c , the contribution of 
charm production to F2, is shown in Fig. [3l The treatment of heavy flavor 
contribution to F2 is constantly evolving [T7] and even within the same 
scheme there are large differences in the expectations of various groups, 
which must affect the uncertainty on the gluon distribution. 




Fig. 3. Compilation of F| c as a function of Q 2 for different values of x, compared 
to expectations of two particular PDFs obtained with the fixed flavor scheme. 



In the DGLAP evolution equations, large gluon densities are expected 
at low x. If the gluon density per unit area becomes large, the probability of 
recombination of gluons becomes high and nondinear effects in the evolution 
may show up [18] which would slow down the increase of the gluon density. 
In addition coherent effects may be expected [19J. It is quite clear by now, 
that the HERA data do not have the power to rule out or confirm the 
presence of such effects. However, there are indications in the HERA data 
that the low x regime of DIS scattering may not be fully described by the 
DGLAP dynamics. This is exemplified by the large fraction of diffractive- 
like events as well as by the observation of an excess of forward jets in the 
low x regime. 
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4. Diffractive scattering in DIS 

If particle production in hard scattering is driven by the DGLAP evo- 
lution, large rapidity gaps (LRG) between hadrons are exponentially sup- 
pressed. Yet, at HERA, about 10% of the DIS events have a large rapidity 
gap, which separates the initial proton from the rest of the hadronic final 
state [20] . It is conceivable that the origin of these LRG is soft in nature and 
that their presence is accounted for in the initial conditions of the evolution, 
which would preserve the validity of the DGLAP dynamics in sufficiently 
hard processes. 

The presence of LRG at high energy is associated with diffractive scat- 
tering which, in soft hadron-hadron collisions, is believed to be due to the 
exchange of a Regge trajectory [2TJ with vacuum quantum numbers - the 
Pomeron (iP). The JP trajectory is characterized by two parameters, the 
intercept a P (0) and the slope a' p . Both have been determined in soft in- 
teractions and the most commonly cited values are 0^(0) = 1.08 [22] and 
a' P = 0.250 GeV -2 [23]. The intercept drives the energy dependence of the 
total cross section, <7 t ot oc s a p(°)^ 1 . The value of a P was later updated [24] 
to a F (0) = 1.0961°;°^. 

The measurements of the diffractive cross section in DIS can be used 
to extract the value of a P (0) |251 126j . A compilation based on the ZEUS 
measurements [26] is presented in Fig. |H The extracted value of cvjp(O) 
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Fig. 4. The value of a w (0) extracted from diffractive DIS scattering as a function 
of Q 2 .The dashed line with the yellow band represents the average value with its 
uncertainty. 

is shown as a function of Q 2 and no dependence on Q 2 is observed. The 
averaged value of a P (0) in DIS tends to lie above the one derived from soft 
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hadron-hadron interactions. A larger value of a P (0) would imply a faster 
growth of the cross section with energy and diffractive scattering in DIS a 
harder process than in hadron-hadron scattering. 



5. QCD factorization in diffraction 

The diffractive contribution to the inclusive structure function Fi is 
called F® . In addition to the usual DIS variables, also depends on 
variables which describe the diffractive final state, 

t = (p-p')\ (2) 

x P = q -^l, (3) 
q-p 

8 = = — . (4) 

P 2q-(p-p') x P {> 

Here p and p' respectively denote the proton four-vector before and after 
the scattering, q is the four-vector of the virtual photon 7* with Q 2 = 
—q 2 . The variable t is the square of the four- momentum exchanged in the 
proton vertex. The variable x P is the fractional proton momentum which 
participates in the interaction with 7* (sometimes denoted by £), 8 is the 
equivalent of Bjorken x but relative to the exchanged object and x P ■ 3 = x. 
QCD factorization is expected to hold for [271 123 EH] and it may then 
be decomposed into diffractive parton distributions (DPDF), f®, in a way 
similar to the inclusive F2, 

dF^(x,Q 2 ,x P ,t) ^ f x P dfP(z,(i,x P ,t) ~ x 2 

-J^u -2^ L dz Jx~~lt ^2,i(-,y (5) 

LLJU P til . J (J yJjJu P iXb £ 

where Fn is the universal structure function for DIS on parton i, \i is 
the factorization scale at which fP are probed and z is the fraction of 
momentum of the proton carried by the diffractive parton i. Diffractive 
partons are to be understood as those partons in the proton from which the 
scattering leads to a diffractive final state. The DGLAP evolution equation 
applies in the same way as for the inclusive case. For a fixed value of x P , 
the evolution in x and Q 2 is equivalent to the evolution in 3 and Q 2 . 

If, following Ingelman and Schlein [30], one further assumes the validity 
of Regge factorization, F® may be decomposed into a universal IP flux and 
the structure function of the IP, 

^^f^ = fjP /p (x P ,t)Fr((3 lQ 2 ) , (6) 
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where the normalization of either of the two components is arbitrary. It 
implies that the x P and t dependence of the diffractive cross section is 
universal, independent of Q 2 and 0, and given by 

2ajp(0)-l 



fjp/p{x P ,t) ~ — 



1 



,(6^-2a^ In x P )t 



(7) 



one of the expectations which agrees with the measurements as discussed 
in the previous section. 

In this approach, the mechanism for producing LRG is assumed to be 
present at some scale and the evolution formalism allows to probe the un- 
derlying partonic structure. The latter depends on the coupling of quarks 
and gluons to the Pomeron. 

An example of diffractive parton distributions derived by the HI exper- 
iment is shown in Fig. [5] [25J. As in the case of the proton PDF, the gluons 




H1 2006 DPDF Fit A 
□ (exp. error) 
EZi (exp.+theor. error) 



0.2 0.4 0.6 0.8 



H1 2006 DPDF Fit B 
- (exp.+theor. error) 



Fig. 5. The singlet, E, and the gluon distribution in diffractive parton distributions 
as a function of the fraction z of the momentum of the colorless exchange responsi- 
ble for the LRG and for different Q 2 values. Two different fit results are presented 
as described in the figure. 

are not well constrained by the inclusive measurements. They are further 
constrained by diffractive jet [31j and heavy flavor [32J[33] production. The 
fact that all diffractive DIS data can be accommodated by the same set of 
DPDF lands support to QCD factorization. 

Gluons constitute a large fraction (about 70%) of the diffractive ex- 
change, implying that DIS induced by gluons is more likely to lead to a 
diffractive final state. In the black-body limit the fraction of diffractive 
events cannot exceed 50% [341 135j. At HERA, this limit is far from being 
reached, however the dominance of gluons in diffractive scattering indicates 
that the onset of unitarity effects may first show up in the gluon. 



ha-text printed on June 3, 2009 



9 



6. QCD factorization breaking in diffractive scattering 

QCD factorization in diffractive scattering is expected to break down 
in hadron-hadron interactions [27] . Indeed, the measurements of diffractive 
dijet production in pp interactions at the Tevatron [36] indicate that their 
rate is by factor 5 to 10 lower than expected from the DPDF extracted 
at HERA. This is shown in Fig. [6j where the relative rate of diffractively 
produced dijets is shown as a function of /3. The suppression of hard diffrac- 




Fig. 6. The ratio of diffractively produced dijets to inclusive dijets as a function of 
j3 in pp collisions compared to expectations of various DPDFs, as indicated in the 
figure, extracted from the HERA data. 

five scattering in pp is understood as due to rescattering of partons which 
did not directly participate in the diffractive process. As a result of the 
rescattering, the LRG is destroyed. This is depicted schematically in Fig. [71 
Some suppression of LRG events is also expected in photoproduction (7p) 




Fig. 7. Left, diagram for diffractive dijet production in pp scattering. Right, dia- 
gram for diffractive dijet production with rescattering. 

at HERA [37], especially in the regime of resolved photon contributions. In 
that respect, the data is not conclusive as the HI experiment observes a 
large suppression for both the resolved photon (more hadrondike) and the 
direct photon (more DISdike) contributions [38] while the ZEUS experiment 
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sees very little suppression if at all [39] . The explanation may lie in the pre- 
liminary measurements of the HI experiment [30] where the cross section 
for diffractive dijet production in *yp is compared to expectations (assuming 
QCD factorization) as a function of the highest transverse energy of the two 
jets, Ej*" 1 . As shown in Fig. [8l the suppression of LRG events is strongest 
at low E^ 1 . This is an interesting observation which indicates that the se- 



H1 PRELIMINARY 

H 1 HERA 99-00 e+ Data f N LO-FR » (1 +S ) 




Ef (GeV> 

Fig. 8. The ratio of dijet diffractive cross section to NLO expectations obtained 
with various DPDFs assuming the validity of QCD factorization as a function of 
the highest transverse energy of the two jets, Ej! . 

lection of higher Et jets renders rescattering less probable. In fact, a closer 
inspection of the CDF results (see Fig. [6]) leads to a similar conclusion. The 
LRG suppression seems stronger at high (3, that is lower diffractive masses 
and therefore lower Et jets. The selection of higher Et jets may be a way of 
selecting a (transversely) smaller partonic configuration of the dissociating 
particle, for which by virtue of color screening the probability of rescattering 
becomes smaller. 

The comparison of the diffractive measurements at HERA and at the 
Tevatron leads to the following puzzle. Of the order of 10% of events ob- 
served in ep collisions cannot be accounted for in pp collisions, yet inclusive 
factorization seems to be preserved (universality of proton PDF). The only 
explanation is that the mechanism of rescattering somehow preserves the to- 
tal cross section. This must have implications for multi-parton interactions, 
one of the big unknowns at LHC. 

7. Hard diffraction and the 3D structure of the proton 

A small fraction of diffractive events in ep scattering at HERA consist of 
exclusive vector meson or single photon (deeply virtual Compton scattering 
DVCS) production. At the high center of mass energies of HERA and 
in the presence of a large scale, these exclusive processes are believed to be 
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mediated by two gluon exchange |41} H2] . The cross section for the exclusive 
processes is expected to rise with center of mass energy W, with the rate 
of growth increasing with the value of the hard scale. This is illustrated in 
Fig. [9] where the W dependence of the 7p cross section for exclusive p, (f>, 
J/ip, ip(2S) and T production is shown together with the W dependence of 
the total 7p cross section (for a detailed discussion and references see f53]). 
The W dependence of the latter is typical of soft interactions. Also shown in 
the figure are the values of the logarithmic derivative 5 obtained by fitting a 
W s dependence. Starting from exclusive J /if) production there is a definite 
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Fig. 9. Total and exclusive vector meson photoproduction cross sections, as a func- 
tion of W. The curves are fits of the form ~ W s 

change in the W dependence, with the tendency for 5 to become larger as 
the mass of the vector meson increases. This is interpreted as the evidence 
of gluon participation in the interaction, since the higher the scale at which 
the gluons are probed the faster the rise with W. 

A further compilation of logarithmic derivatives [13] , which also includes 
the measurements of exclusive processes in DIS, for p, <f>, J/if> and DVCS 
as a function of scale defined as Q 2 + My, where My is the mass of the 
vector meson, is presented in figure [TU1 With increasing hardness of the 
interaction, the t distribution is expected to become universal, independent 
of the scale and of the final state. This is because the harder the scale 
the more pointdike the virtual photon becomes. The exponential slope of 
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Fig. 10. Left: logarithmic derivatives S = dlog o~("f* p) / dlog W as a function ofQ 2 + 
M 2 for exclusive processes. Right: Exponential slope of the t distribution measured 
for exclusive processes as a function of Q 2 + M 2 . 



the t distribution, b, reflects then the size of the proton. A compilation of 
measured b values (see [44] for details and references therein) is presented in 
figure [TU1 Around Q 2 +M 2 of about 15 GeV 2 the b values becomes universal. 

Since exclusive vector meson production or DVCS at HERA energies 
is driven by gluons, the b value probes the size of the gluon cloud in the 
proton [45] . The universal value of b ~ 5 GeV -2 corresponds roughly to 
a radius of 0.6 fm which is smaller than the 'charge' radius of the proton 
which is 0.8 fm. This is the first indication that the gluons are well con- 
tained within the proton. This finding has implications for multi-parton 
interactions. In particular the more central the pp interaction the higher 
the probability of a second partonic interaction [45] . 



8. Hadronic final states 

At LHC, much of the background estimates rely on proper simulations 
of the hadronic final states, of which multi-parton interactions are only one 
of the elements. 

Most of the MC generators base the formation of the hadronic final sate 
on the properties of the DGLAP dynamics, for which subsequent emissions 
of partons, before and after the hard interaction, are strongly ordered in 
transverse momentum kx, with &t increasing steadily while the momentum 
fraction x decreases. At low x, where effects of the BFKL [36] dynamics 
may be felt, the ordering in )zt may be lost while subsequent emissions 
are strongly ordered in x. The two different mechanisms are schematically 
depicted for DIS in Fig. [TT1 Note that the contribution of large higher 
order corrections to the DGLAP evolution may also appear as loss of )zt 
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ordering. The loss of strong kx ordering may result in the appearance of 
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Fig. 11. Schematic representation of the DGLAP (left) and BFKL (right) parton 
radiation patterns 



high Et forward (in the proton direction) jets, of jets with > Q and a 
decorrelation in the azimuthal angle between jets. All of these effects have 
been observed at HERA. 

The cross section for producing forward jets, defined by 1.7 < rf 01 < 2.8, 
Xjet > 0.035 and with transverse momentum pr fulfilling 0.5 <p\IQ 2 < 5, 



as measured by the HI experiment [17], is shown in Fig. [12] as a function 
of x. The measurements are compared to various MC generators. A clear 
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Fig. 12. Cross section for producing forward jets, with 1.7 < ?f Qt < 2.8, Xj Ct > 0.035 
and with transverse momentum fulfilling 0.5 < p\jQ 2 < 5 as a function of 
Bjorken x, xej- Also shown are expectations of NLO calculations and various MC 
generators. 



excess of forward jets is observed at low x, which cannot be accounted for 
neither by higher order DGLAP corrections nor by any other MC genera- 
tor. The highest rate is expected in the Color Dipole Model (CDM) [38] 
generator in which the lack of )zt ordering is implicit in the way the color 
dipoles emit radiation. If the DGLAP dynamics, as represented by the RG 
generator [49J , is supplemented by the addition of a resolved virtual photon 
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Fig. 13. Schematic representation of diagrams with direct (DIS) and resolved 7* 
contributions included in the simulation of DIS hadronic final states. 

contribution [50J, the rate of forward jet production increases but still fails 
to describe the lowest x measurements. The addition of a resolved virtual 
photon contribution is supposed to mimic the lack of kx ordering as shown 
schematically in Fig. [T3l 




Fig. 14. Cross section for dijet production as a function of x J 7 0ts , the estimator of 
the 7* momentum carried by the two jets, for various ranges of transverse energy 
measured in the 7*p center of mass system, Ej>, and Q 2 . Also shown are predictions 
of NLO calculations with and without the resolved 7* contribution. 

The addition of a resolved photon contribution also helps in describing 
the cross section for dijet production in DIS, in particular in the regime 
where the E 2 -, of the jets is higher than Q 2 |51j . This is shown in Fig. [TJ] 
where the dijet differential cross section is plotted as a function of x^ ts , the 
estimator of the 7* momentum carried by the two jets, for various ranges of 
transverse energy measured in the j*p center of mass system, Ej,, and Q 2 . 
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However even the addition of the resolved photon contribution falls short of 
describing the low Q 2 region where x is lowest. 

The azimuthal correlations in multijet production have been measured 
both by the HI [52J and the ZEUS [53J experiments. The NLO calcula- 
tions of dijet production fail to describe the rate of jets with the difference 
in azimuthal angle A0 < 120° as shown in Fig. [T5l in particular at low 
x. However, an addition of a 3 terms seems to cure the problem and also 
properly describes the rates for trijet events. 




Fig. 15. DIS cross section for dijet and trijet production, when the two highest Et 
jets are separated by an azimuthal angle A</> < 120°, as a function of Bjorken x - 
x-By Also shown are the expectations of NLO calculations without and with 0(a^) 
corrections. 

All these findings strongly suggest that at low x the description of 
hadronic final states requires corrections even beyond NLO DGLAP. 



9. Summary 

The large fraction of DIS diffractive events, mainly driven by gluons, 
could suggest the approach to the unitarity limit, in which case non-linear 
effects could bias the measurements. The observation of an excess of for- 
ward jets in the low x regime would suggest that the transverse momentum 
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ordering of gluon radiation, inherent in the DGLAP evolution, breaks down 
at low x. The same conclusion may be derived from the appearance of the 
so called 'resolved photon' contribution necessary to described DIS dijet 
production in the low x regime. At HERA, all these effects constitute a rel- 
atively small fraction of the total cross section (typically 5%) and therefore 
may not be visible in the fully inclusive measurements. At LHC, these low x 
effects may be greatly enhanced. In particular, the increased jet multiplicity 
and the shift of the high transverse momentum activity toward the proton 
direction could result in events with missing transverse energy greater than 
expected from Monte Carlo generators. 

The non-uniform space distribution of partons in the proton, partly 
covered by the newest generation of MC programs, and the strong break- 
ing of diffractive QCD factorization due to rescattering, with potentially 
non-trivial dependence on kinematics, may lead to surprises in modeling 
the underlying event and multi-parton interactions accompanying hard pp 
interactions at the LHC. 
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